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Though QCD Factorization(QCDF) could produce sufficiently large B — > Kn branching ratios 
close to experiments, the predicted direct CP asymmetry in i?" — > K^-k~ decay is however still quite 
below experiment, even with a large negative phase in the annihilation term, according to existing 
calculations. This suggests the presence of an additional strong phase in the decay amplitude which 
could come from the long-distance final state interactions(FSI) like charming penguin. In this paper, 
we show that, by adding to the B Ktv decay QCDF amplitude, a real part and an absorptive part 
with a strength 10% and 30% of the penguin amplitude, respectively, we could bring the B —> Kn 
branching ratios and the — > K~-k^ CP asymmetry close to the measured values. We also find 
that the color-allowed electroweak penguin is appreciable and that from the QCDF electroweak- 
strong penguin interference terms and the measured B~ — > K^tt' we obtain (9.0 ± 0.3) x 10~® 
for the i3° K^tv'^ branching ratio, a bit lower than experiment. Similarly, the predicted value 
for B~ if^Tr" in terms of the — > K~tv'^ branching ratio agrees well with experiment. This 
suggests that the measured value for — > K %o in fact should be lower than the current measured 
value. 

PACS numbers: 13.25Hw, 12.38.Bx 



in: 

> 
O 

m 

00 

o 
o\ 
o 



X 



I. INTRODUCTION 



Existing calculations of the penguin-dominated B — > 
Kit decays in QCD Factorization(QCDF) [J, [2, [1| seem 
to produce B Kit branching ratios mre or less in agree- 
ment with experiments, with a moderate annihilation 
contribution, but the predicted direct CP asymmetry in 
B^ — !■ K^TT^ decay is smaller than experiment by a fac- 
tor of 2 or more, even with a large negative strong phase 
for the annihilation terms as in scenario S4 [3| . In fact in 
QCDF, and also in naive Factorization model ^4,], the ab- 
sorptive part(the imaginary part) of the short-distance 
Wilson coefficients in the non leptonic decay effective 
Hamiltonian are negative and of the same sign as the 
real part, the direct CP asymmetry for B^ tt+tt" and 
B° if ~7r+ are negative and positive, respectively, in 
opposite sign with measurements [1]. To reverse the sign 
of the CP asymmetry and to produce a large CP asymme- 
try, one would need a large annihilation term with a large 
negative phase. Another possibility is a long-distance fi- 
nal state interaction(FSI) term, probably mainly absorp- 
tive, to generate a large strong phase, but not to overes- 
timate the branching ratios, like charming penguin due 
to B ^ D,D ^ Ktt and B ^ D*D* Ktt inelastic 
FSI effect d 0, B i, E [H which are CKM-favored, 
color-allowed processes and could produce a large ab- 
sorptive part in i? ^ Kit decays. In this paper, we 
show that by adding to the QCDF amplitude, a penguin- 
like amplitude with a real and absorptive part with a 
strength 10% and 30%, respectively, of the penguin am- 
plitude, we could obtain the B —^ Kt: branching ratios 
and the B^ — > K~Tr^ CP asymmetry close to experi- 
ments. We also find that the color-allowed electroweak 
penguin(EW) from the operators O7 and O9 contribute 
appreciably through the strong penguin- EW penguin in- 



terference terms in B^ K^tt^ and B^ A'°7r° de- 
cays. Then using only the QCDF expressions for the 
color-allowed tree-penguin(T-P), color-suppressed tree- 
penguin(C-P), color-allowed electroweak-penguin(EW- 
P) interference terms and the measured B{B~ — > K^tt^) 
and S(B° ^ K-TT+), we could predict B{B- K-jt^) 
and B{B° K'^tt^). We find that the predicted B{B~ -> 
K^tt'^) agrees well with experiment while the predicted 
B{B'^ — > T^^tt") is below the measured value. This 
strongly suggests that the measured B{B'^ — > X^tt") is a 
bit too high as shown in the following. In the next section 
we will present numerical results for the B — )■ Kt: branch- 
ing ratios and CP asymmetry in QCDF factorization with 
the addition of a small penguin-like amplitude to gener- 
ate a CP asymmetry for _B" K^tt^ decay. We then 
give prediction for B{B- R-jr") and 6(5" -> K^tt^) 
from the electroweak penguin and strong penguin inter- 
ference terms. 



II. B -* K-K DECAY IN QCD FACTORIZATION 

The B M1M2 decay amplitude in QCDF is given 

by 

A{B MM - ^ E ^P>^Ks X 

^ p—u,c 

/ 10 10 \ 

- ar(MiM2|0.|B)H + E Ib/mJalA , (1) 



where the QCD coefficients contain the vertex 
corrections, penguin corrections, and hard spectator 
scattering contributions, the hadronic matrix elements 
{MiM2\Oi\B)H of the tree and penguin operators Oi are 
given by factorization model [J, Il2l | , bi are annihilation 
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contributions. The values for a^-,p = u^c , computed 
from the expressions in [l], Q at the renormahzation scale 
fi — rrif,, with mt, = 4.2 GeV, are : 

a1 = -0.033 - 0.013 i + 0.0009 pH exp(i0_f/), 
= -0.027 - 0.017i + 0.0009 pH exp{i(l)„), 



a^e = -0.045- 
= -0.0004 



0.003 1, < 
- O.OOOH, 



-0.042 - 
= 0.0004 



O.OlSi, 
- 0.0001 i, 



no 



= -0.0011 - 0.0001 i - 0.0006 exp(i0_f/) 



a^'o = -0.0011 



0.0006 i- 0.0006 pHexp{i(j)„). (2) 



u _ „P _ 



for i ~ A,6, 8, 10. For other coefhcients, a 

ai ^ 1.02 + 0.0l5i - 0.012 PH exp{i(j)H), 
02 = 0.156 - 0.089 i + 0.074 pjjexp(ii<^//), 
ag = 0.0025 + 0.0030 i - 0.0024 pH exp(i0jj ), 
as = -0.0016 - 0.0034 i + 0.0029 pH exp{i(l)H), 
aj = -0.00003 - 0.00004 i - 0.00003 pn exp(i0_f/) 
ag = -0.009 - 0.0001 i + 0.0001 pn exp(i0_ff ). (3) 

where the complex parameter pn exp(i(f)H) represents the 
end-point singularity contribution in the hard-scattering 
corrections Xh = {1 + pn exp{i(t>H)) Ht^) [EH- 
For the annihilation terms, we have : 

62 = -0.0041 -0.0071pAexp{i(j)A) - 0.0019{pA exp(^0A))^ 

63 = -0.0071 -0.016pA exp{i(j)A) - 0. 0093{p a exp{i(P a))^ , 
6^'" = -0.00012- 0.00016 PA exp(i0A) 

-hO.OOOOO3(pAexp(i0A))2. (4) 



where bi are evaluated with the factor fBfMifM2 
included and normalized relative to the factor 
fKF^^{m\ — m^) in the factorisable terms, and 
Pa , like pu, appears in the divergent annihilation term 

XA^{l + PAexp{ld,A))H7^)- 

The B Ktt decay amplitude with the factorisable 
part 4] and the annihilation term [l], [3, [13] are: 

AiB- ^ K-n') = -^^f,,F,^^ml){m% - ml) 
{VubV;,ai + {VutV;, + VcbV;,)[ai + + {ae + as)r^]) 



B 



X [VubV:,a2 + (KbK*. + VcbV;,) X -(ag - aj) 
.Gf f , f 

-^ — JBJKJtt 

X [VubV:M + (KbK. + V,bV:,) X (63 + 6™)] 



(5) 



A{B- 



X [Vubv:,h + {Vubv:, + Vcbv:,) x (63 + )] 



(6) 



and for : 



A{B° 



G 



KblCsfli +(K6K*s +K6K*s)[«4 +aio +(a6 + aB)r^ 
Gp 



{Vubv:, + v,bv:,) X (63 



"3 



(7) 



A{B' ^ K\') = ^^f^F,^^ml){ml - ml] 



xiVubV:, + v,bV;,) 

. G F 



1 .In 



^—UFi^{ml){ml~ml) 



VubV:,a2 + {VubV:, + KfcKl) X -(ag - a,) 



+i^fBfKfw 



where 



{Vubv:, + v,bv:,) X (63 



"3 • 



(8) 



r-v — 7 tt'' — : 7 is the chirally-enhanced 

terms in the penguin Og matrix element. We also need 
the B'^ TT^TT^ amplitude: 

A{B- ^ n-n") = -^^^^^„^'^(m^)(m| - ml) 



X -(ag - ay + aio + asr^) 



(9) 



We see that the B Ktt decay amplitudes consist of 
a QCD pcnguin(P) 04 + aer^ , a color-allowed tree(T) 
ai, a color-suppressed tree(C) a2 , a color-allowed elec- 
troweak penguin(EW) ag — ay, a color-suppressed elec- 
troweak penguin(EWC) aio -t- asr^ terms. Because of 
the relative large Wilson coefhcients, the QCD penguin, 
the color-allowed tree and the color-allowed elcctrowcak 
contributions give the major contribution in B Kn 
decays. 

For the CKM matrix elements, since the inclusive and 
exclusive data on \Vub\ differ by a large amount and 
the higher inclusive data exceeds the unitarity limit for 
Rb = \VudV*b\/\VcdV*b\ with the current value sin(2/3) = 
0.681±0.025 1J|, we shall determine \Vub\ from the more 
precise \Vcb\ data. We have p^ : 



|K.H^^|sin/?WH- 



IK 



ud I 



• 2 

sm a 



(10) 



With a = (99lg^)° [li] and \Vcb\ = (41.78±0.30±0.08) x 
10-3 [ll], we find 



Kb I = 3.60 X 10" 



(11) 



in good agreement with the exclusive data in the range 
I Kb I = 3.33 - 3.51 • A recent UT fit also gives 
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IKftl = (3.60 ±0.12) X 10-3 [T3|. The measurements 
of the Bs — Bg mixing also allow the extraction of \Vtd\ 
from Bd — Bd mixing data. The current determination 
[H] gives \VtdlVts\ = (0.208t5]:J]!]^) which in turn can be 
used to determined the angle 7 from the unitarity relation 



\Vtd\ = „,,f' |sm7 



tb\ 



1 + 



• 2 

sm a 



(12) 



with \Vtb\ = 1, we find 7 = 66° which implies an angle 
a — 91.8°, in good agreement with the value found in the 
current UT-fit value of (88 ±16)° In the following 
in our B Ktt decay calculations, we shall use the uni- 
tarity triangle values for \Vub\ and 7. For other hadronic 
parameters we use the values in Table 1 of '5] and take 
TO,(2GeV) = SOMeV. For the B ^ tt and B ^ if 
transition form factor, we use the current light-cone sum 
rules central value 20] : 



Fo^^(O) = 0.258, F(f (0) = 0.33 



(13) 



The computed branching ratios and direct CP asym- 
metries, with ph = ^,(t'H — and c/ja = —55° as in 
scenario S4 of [1] are shown in Fig[T] and Fig[2] as func- 
tion of Pa ■ For convenience we also give in Table U and 
Table nil the computed values a,t pA — ^ as in S4 with and 
without the additional penguin-like contribution SP. 



Modes 


SP^O 


SP^O 


Exp [21] 


B- TT-n" 


5.7 


5.7 


5.59 ±0.4 


B- -* K-Tv^ 


10.3 


12.7 


12.9 ±0.6 


B- K°TT- 


18.1 


22.7 


23.1 ± 1.0 


B° ^ K-7V+ 


15.5 


19.7 


19.4 ±0.6 




6.8 


8.8 


9.8 ±0.6 



TABLE I: The CP-averaged B Kn Branching ratios in in 
unit of 10~^ in QCDF with and without additional penguin- 
like contribution SP and with pA = 1.0, (j)A ~ —55° 



Modes 


5P = 


SP^O 


Exp [21] 


B- TT-yr" 


0.0 


0.0 


0.06 ±0.05 


B' K'Tv'^ 


0.02 


-0.06 


0.05 ± 0.025 


B- K°n- 


0.004 


0.01 


-0.009 ± 0.025 


BO ^ K-7V+ 


-0.01 


-0.08 


-0.098^H1^ 


50 ^ K%'' 


-0.02 


0.01 


-0.01 ±0.0.10 



TABLE II: The direct B ^ Ktt CP asymmetries in QCDF 
with and without additional penguin-like contribution SP and 
with PA = 1.0, 4>A — —55° 



is much smaller than the measured values. Considering 
various theoretical uncertainties, one could say that the 
real part of the B — > Kn decay amplitudes is more or less 
well described by QCDF, but the absorptive part of the 
amplitude needs some additional contribution to account 
for the observed CP asymmetries. As mentioned earlier, 
the factorisable contribution in the QCDF amplitude give 
the wrong sign for the CP asymmetries in B'^ — s- K~Tr~^, 
one needs a sizeable absorptive part to reverse the sign of 
the asymmetries. As shown in Fig[T]and FigO we would 
need an annihilation term with a very large pA to pro- 
duce correct B Ktt branching ratios and i?" K^tt'^ 
asymmetry. This leaves us with charming penguin, as 
mentioned earlier, as a possible contribution to the large 
strong phase in i? Ktt decays. This possibility is con- 
sidered in this paper by adding to the QCD penguin a 
small term: 



SP — (ra4 ± raQr^){di ± id2) 



(14) 



where ra4 and rag are the real part of 04 and ag, respec- 
tively. We note a previous fit with charming penguin 
for B — > PP, PV in . The results shown above are 
obtained with di ~ 0.08 , ^2 — —1/3. We see that the 
branching ratios for B'^ K^tt^, B^ — > K^tt'^ and 
B~ K°TT- and the CP asymmetry for B° K-tt+ 
are in good agreement with experiment. The B{B'^ —> 
K'^tt'^) is slightly below the measured value. Our results 
seem reasonable, as the predicted B{B^ — > tt~tt'^) which 
is sensitive F{^^{Q) agrees well with the measured value. 

Thus QCDF with mainly absorptive additional pen- 
guin contribution seems to describe rather well the B 
Ktt observed branching ratios and the B^ A'^7r+ 
direct CP asymmetry. As all the B — > Ktt branch- 
ing ratios have been measured with an accuracy at the 
level of 0(10"^), one could look at the difference in 
branching ratios to test QCDF predictions for the elec- 
troweak and tree- level contributions. The color- favored 
electroweak penguin is CKM-favored, so its interference 
with the strong-penguin terms is about the same size 
as that of the color-favored tree and penguin interfer- 
ence. For example, from Eq.(l5][8]), one sees that the 
differences _2r(S- ^ K-tt") - r{B° K'tt+) and 
2T{B° K"tt°) - T{B- -> K"tt-) are essentially given 
by the color-allowed electroweak and strong penguin in- 
terference term, so with QCDF expression for the inter- 
ference term, one could predict B{B''^ — > A'°7r°) in terms 
of B{B~ K^TT~) or vice versa. We have: 



B{B° ^ K\") = [B{B- - 
B{B- K-TT°) = [6(5° 
B{B° K-TT+) = [B{B- 



K\-)rb + SBi]/2 (15) 
> K-TT+) + SB2]/2rh (16) 
^K\-)n + 6B3] (17) 



As can be seen, without the additional penguin contri- 
bution, the computed branching ratios are slightly below 
experiments, but the CP asymmetries for — > K^tt^ 



and where SBi, SB2, SB3 are respectively the com- 
puted difference 2B{B° K°tt°) - nB[B- K°tt-), 
2nB{B- -> K-tt") - B{B° K~TT+), B{B° 
K-TT+) - rbB{B- -> K^TT-) and n = tbo/tb-- We 
find, in unit of 10^®, with pA = 1,0A = —55° as with 
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FIG. 1: The computed and measured CP average branching ratios. The horizontal line are the measured values [2l|] with the 
gray areas represent the experimental errors, (a), (b), (c), (d) in the left and right figure represent the values for B~ K~-k^ , 
B- K°7v-, B° K--K+ and B^ ^ K°-k° respectively. The curves (al)-(dl) and (a2)-(d2) are the corresponding QCDF 
predicted values for 4>a ~ —55°, without and with additional penguin contribution respectively. 




FIG. 2: The same as in Fig[l]but for the computed CP asymmetries. 



the computed values in Table [J and Table [TTlf only leading 
terms are shown): 

5Bi = -3.52 [-3.84(EW) + 1.42 cos(7)(C)] (18) 

5B2 = 3.97 [3.84(EW) - 1.42cos(7)(C)] (19) 
5B3 = -1.52 [-5.53cos(7)(TP) 

-0.65cos(7)(TA) + 0.69(EWC) + 0.36(TT)] (20) 

In the above expressions, the main contributions are 
from the interference of the color-allowed tree, color- 
suppressed tree, and electroweak penguin with the strong 
penguin terms and are shown inside the brackets in the 
above expressions. Then from Eq. (|15m7| . we obtain (in 
unit of 10"*^) 

S(B° ^ /fV) = 9.0±0.3, (21) 
S(B- ^ if"7r°) = 12.5±0.3. (22) 
^(5° ^ /f"7r+) = 20.0 ±0.6, (23) 

which agrees with experiment, to within the current 
accuracy. In particular the agreement with experiment 



for B{B K TT^) shows that the electroweak-penguin 
interference observed in _B — > Kir decay rates difference is 
well described by QCDF. The predicted B{B'^ K^tt°) 
is a bit lower than experiment, but this prediction is 
on firm ground since dBi « —dB2 and the predicted 
B{B^ K~Tr°) in terms of dB2 agrees well with exper- 
iment, as mentioned above. In fact the discrepancy be- 
tween experiment and the above prediction for B{B'^ —> 
K'^tt'^) could be understood in the following way : First 
of all, Eqs. (|18ll^ shows that the electroweak-penguin as 
well as the color-suppressed tree with the strong penguin 
interference terms in B{B^ K^tt'^) and in S(i?° 
K'^tt^) are essentially the same in magnitude but opposite 
in sign. This can be seen from the decay amplitudes given 
in Eqs.([5]-[5]). Then these interference terms cancel out 
leaving only the color-allowed tree-penguin interference 
term in the sum of r{B- K'tt^) and T{B" K°n°). 
Since the color-allowed tree-penguin interference terms, 
the penguin and annihilation terms are essentially the 
same in the sum 2{r{B- R-n") + T{B" K°7r°) 
and {r{B- K^n- + T{B° K-tt+)), by neglecting 
small quadratic term in tree and electroweak penguin 
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contributions, we obtain the approximate relation: 

2[B{B- ^ K-nyt + ^ 

= [B{B- K°TT-)rb + B{B° K-tt+)] (24) 

This relation has been first given in [2^ [l^j, but was 
derived independently in [25l | from the isospin B —> Ktt 
amplitudes: 



A 



{r{B- i^°7r-) + r(B° - 

2 [r{B- -> K-7r°)+T{B" 
4 
3 



^|i?3p - ^Re(i?*Bie'^) 



V3 



K-7r+) 
(Crso) (25) 



where Bi^B^ are respectively the isospin / = 1/2 and 
1 = 3/2 tree-level B Ktt amplitudes. With the mea- 
sured branching ratios, we get (43.69 ± 2.4) x 10~^ and 
(40.9 ± 2.0) X 10"^ for the l.h.s and r.h.s of Eq.dM]), and 
a difference 2.72 x 10~^ between the two values. Ignor- 
ing the errors, we see that this difference is rather large, 
of the size of the tree-penguin interference term, while, 
theoretically, it should be of O(10~^) or less and has a 
value given by dBi+SB2 — 0.43 x 10~^. The problem, we 
suspect, seems to be the measured S(i?° — + K'-'tt^) which 
has been lowered over the years. It was (11.5±1.0) x 10~^ 
di] and exceeds the upper limit of (10.8 ± 0.3) x 10"*^ 
in 2005, as shown in Eq.(IT5| (the old Ktt puzzle [13]). 
With our value for S(i?° — ^ K'^tt'^) , the difference would 
be reduced to 0.9 x 10~^ and becomes consistent with 
Ea. ([M)) . within the present accuracy. We thus expect 
the experimental value for B{B^ K^tt'^) would be 



lower than the current value and is close to our pre- 
dicted value 9.0 ± 0.3 x 10"^. (The Belle value is in 
fact (9.2±0.7lo ?) x 10"^ HH). 



III. CONCLUSION 

By adding a mainly absorptive, penguin-like contri- 
bution to the B Ktt decay QCDF amplitude, we 
show that QCDF could successfully predict the B Ktt 
branching ratios. We obtain the correct magnitude and 
sign for the B^ K~tt^ CP asymmetry. The B~ 
K^TT^ CP asymmetry is predicted to be of the same sign 
and magnitude as that of B^ K^tt^ which should be 
checked against more precise data. We also show that 
QCDF describes well the electroweak penguin contribu- 
tion in the B~ — > K~tt^ and B^ — > K^tt^ branching ra- 
tios and show that the B^ — > K^tt^ measured branching 
ratio should be lower than the current value. 
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